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Abstract 

In this paper we considcr thc minimum timc population transfcr problem for the z-component 
of the spin of a (spin 1/2) particle driven by a magnètic field, controlled along the x axis, with 
boundcd amplitudc. On the Bloch sphere (i.e. after a suitable Hopf projection), this problem can 
be attacked with techniques of optimal syntheses on 2-D manifolds. 

Lct (—E,E) be the two energy levels, and |0(í)| < M the bound on the field amplitudc. For 
each couplc of vàlues E and M . we determino thc timc optimal synthcsis starting from the levcl —E 
and we provide the explícit expression of the time optimal trajectòries steering the state one to the 
state two, in terms of a parameter that can bc computed solving numcrically a suitable equation. 

For M/E << 1, every time optimal trajectory is bang-bang and in particular the corresponding 
control is periòdic with frequeney of the order of the resonance frequeney u>r = 2E. 

On the other side, for M/E > 1, the time optimal trajectory steering the state one to the state 
two is bang-bang with exactly one switching. Fixed E we also prové that for M — > oo the time 
needed to reach the state two tends to zero. In the case M/E > 1 therc arc timc optimal trajectòries 
containing a singular arc. 

Finally we compare these results with some known results of Khaneja, Brockett and Glascr and 
with thosc obtained by controlling thc magnètic field both on thc x and y directions (or with one 
external field, but in the rotating wave approximation) . 

As byproduct we prové that the qualitative shapc of the time optimal synthcsis presents different 
patterns, that cyclically alternate as M/E — > 0, giving a partial proof of a conjecture formulated in 
a previous paper. 
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1 Introduction 



1.1 Preliminaries 

The issue of designing an cfficicnt transfer of population between different atòmic or molecular levels is cru- 
cial in atòmic and molecular physics (see e.g. [7, 22, 24, 34]). In the experiments, excitation and ionization 
are often induced by means of a sequence of làser pulses. The transfer should be as efficient as possible in 
order to minimize the effects of relaxation or decoherence that are always present. In the recent past years, 
people started to approach the design of làser pulses by using Geomètric Control Techniques (see for instance 
[13, 19, 20, 27, 31]). Finite dimensional closed quantum systems are in fact left (or right) invariant control 
systems on SU(n), or on the corresponding Hilbert sphere S 2 "" 1 C C", where n is the number of atòmic or 
molecular levels. For these kinds of systems very powerful techniques were developed both for what concerns 
controllability [21, 23, 26, 32] and optimal control [3, 12, 25]. 

The dynamics of a n-level quantum system is governed by the time dependent Schròdingcr cquation (in a 
system of units such that h = 1), 

m 

ix(t) = (H + J2 iïj(t) H j)x(t) (1) 
i=i 

where x(.), defined on [0, T] is a function taking vàlues on the state space which is SU{n) (if we formulatc 
the problem for time evolution operator) or the sphere 5 2 "" 1 (if we formulate the problem for the wave 
function). The quantity Hq called the drift Hamiltonian is an Hcrmitian matrix, that is natural to assumc 
diagonalizcd, i.c., Hq = diag(Ei, E n ), where E\,...,E n are real numbers representing the energy levels. 
With no loss of gener ality we can assume Y^j=i Ej = 0- The real valued controls Oi(.), ...,Cl m (.), represent 
the external pulsed field, while the matrices Hj (j = 1, ...,m) are Hermitian matrices describing the coupling 
between the external fields and the system. The time dependent Hamiltonian H(t) := Hq + Y^JLi Qj(t)Hj is 
called the controllcd Hamiltonian . 

The first problem that usually one would like to solve is the controllability problem, i.e. proving that for 
every couple of points in the state space one can find controls steering the system from one point to the other. 
For applications, the most interesting initial and final states are of course the eigenstates of Hq. 

Thanks to the fact that the control system (1) is a left invariant control system on the compact Lie group 
SU(n), this happens if and only if 

Liejziío, iHi, ...iH m } = su(n), (2) 

(see for instance [32]). The problem of finding easily verifyable conditions under which (2) is satisfied has been 
deeply studied in the literature (see for instance [6]). Here we just recali that the condition (2) is genèric in 
the space of Hermitian matrices. 

Once that controllability is proved one would like to steer the system, between two fixed points in the state 
space, in the most cfficicnt way. Typical costs that are interesting to minimize for applications arc: 

pT to 

• Energy transfered by the controls to the system. / V^fí 2 (í) dt, 

• Time of transfer. In this case one can attack two different problems one with boundcd and one with 
unboundcd controls. 

The problem of minimizing time with unbounded controls is now well understood [1, 27]. On the other side, 
the problems of minimizing energy, or time with bounded controls are hopeless in general. Indeed optimal tra- 
jectòries must satisfy a first order necessary condition for optimality called the Pontryagin Maximum Principio 
(in the following PMP, see for instance [3, 25]) that generalizes the Weierstrafi conditions of Calculus of Vari- 
ations to problems with non-holonomic constraints. For each optimal trajectory, the PMP provides a lift to 
the cotangent bundle that is a solution to a suitable pseudo-Hamiltonian system. Hence, the first difficulty 
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comes from thc problcm of integrability of a Hamiltonian system (that generically is not integrable except for 
very special costs). Second, one should manage with some special solutions of the PMP, the so callcd abnormal 
and singular extremals (see for instance [8]). Finally, even if one is able to find all the solutions of the PMP 
(called extremal trajectòries), it remains the problem of selecting, among them, the optimal trajectòries, that 
usually is even a more difficult problem. For that purpose high order necessary conditions for optimality have 
been studied, like Clebsch-Legendre conditions, higher-order maximum principle, envelopes, conjugate points, 
index theory (cf. for instance [3, 12] and references therein). For these reasons, usually, one can hope to find 
a complete solution to an optimal control problem for very special costs, dynamics and in low dimcnsion only 
([4, 12, 17, 33]). 

In [13, 14, 15, 16] a special class of systems, for which the analysis can be pushed much further, was 
studied, namely systems such that the drift term Hq disappear in the interaction picture (by a unitary change 
of coordinates and a change of controls). For these systems thc controllcd Hamiltonian rcads 



H(t) 



( E x /iiíli (í) ••• \ 

A*aflS(í) '■■ '■■ 

■ '• '■ E n _\ /A n _iíï„_i(í) 

V ••• /*„_!«*_!(*) E n ) 



(3) 



Hcre (*) denotes the complex conjugation involution. The controls fli, . . . , íl^-i are complex (they play the 
role of the real controls ííi, . . . , fi m in (1) with m — 2(n — 1)) and fj,j > 0, (j = 1, . . . , n — 1) are real constants 
describing the couplings (intrinsic to the quantum system) that we have restricted to couple only levels j and 
j + 1 by pairs. 

For n = 2 the dynamics (3) describes the evolution of the z component of the spin of a (spin 1/2) particle 
driven by a magnètic field controlled both along the x and y axes, while for n > 2 it represents the first n levels 
of the spectrum of a molecule in the rotating wave approximation (sec for instance [5] ) , and assuming that 
each external fields couples only close levels. The complete solution to the optimal control problem between 
eigenstates of Hq = diag(Ei, . . . , E n ), has been constructed for n = 2 and n = 3, for the minimum time 
problem with boundeí 
(with fixed final time) 



problem with bounded controls (i.e., |íl, | < Mj) and for thc minimum energy problem J Q T X^=i 1^/(01 2 ^ 



Remark 1 For the simplest case n = 2 (studied in [13, 20]), the minimum time problem with bounded control 
and the minimum energy problcm actually coincide. In this casc the controllcd Hamiltonian is 

m) = ( ~t "í? ), i«i < m, (4) 



íT(í) E 

and thc optimal trajectòries, steering the system from the first to the second eigenstate of Hq = diag{—E, E), 
correspond to controls in resonanec with the energy gap 2E, and with maximal amplitude i.e. íl(í) = 
j^ e i[(2E)t+<t>] ^ w j iere ^ g jq^ 2tt[ is an arbitrary phase. The quantity lür = 2E is called the resonance frequeney. 
In this case, the time of transfer T'c is proportional to the inverse of the làser amplitude. More precisely (see 
for instance [13]), T c = n/(2M). 

For n — 3 thc problcm has been studied in [14, 16] and it is much more complicated (in particular when the 
coupling constants fií and /12 arc diffcrcnt). In thc case of minimum time with bounded controls, it requircs 
some nontrivial tcchnical tools of 2-D syntheses theory for distributional systems, that have been developcd in 
[16]. 

For n > 4 the problcm is hopcless, but in [15] it has been proved that the optimal controls steering the 
system from any couple of eigenstates of Hq are in resonance , i.e. they oscillate with a frequeney equal to the 
difference of energy between the levels that the control is coupling. More precisely 

íl,- = A j (t)e i ^+ 1 - E ^ t+ ^\ ] = 1, ...,n - 1 (5) 
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where Aj(.) are real functions describing thc amplitudc of thc extcrnal fields and <pj are arbitrary phascs. 
Actually, this result holds for more general systems, initial and final conditions, and costs (see [15]). 



The problem of minirnizing time with bounded controls or energy is even more difficult if it is not possi- 
ble to eliminate the drift Ho- This happens, for instance, for a system in the form (3) with real controls 
flj(í) = O* (í) , j = 1, n — 1, as we are going to discuss now. (For more details on the elimination of the drift 
see [13, 14, 15].) 

1.2 A spin 1/2 particle in a magnètic field 

In this paper we attack the simplest quantum mcchanical model interesting for applications for which it is 
not possible to eliminate the drift, namely a two-levcl quantum system driven by a real control . This system 
describes the evolution of the z component of the spin of a (spin 1/2) particle driven by a magnètic field 
controlled along the x axis. Equivalently it describes the first 2 levels of a molecule driven by an external 
field without the rotating wave approximation. The dynamics is governed by the time dependent Schròdingcr 
equation (in a system of units such that H = 1): 

i^=H(t)m (6) 

where ip(.) = {M^M^V ■ [0,T\ -> C 2 , £ 2 =1 |^(í)| 2 = 1 (i.e. ip(t) belongs to thc spherc S 3 C C 2 ), and 

HW=ín,f> Q( ' } ), (7) 



íl(t) E 

where E > and the control f2(.), is assumed to be a real function. With the notation of formula (1), the drift 
Hamiltonian is H = ^ ~^ ^ ^ , whilc Hi = ( J J J , and thc controllability condition (2) is satisfied. 

Notice that for a spin 1/2 system, it is equivalent to treat the problem for the wave function or for the time 
evolution operator since S 3 is diffeomorphic to SU (2). The aim is to induce a transition from the first eigenstate 
of Ho (i.e., \tpi | 2 = 1) to any other physical state. We recali that two states ip, ip' £ S 3 are physically equivalent 
if thcy differ for a factor of phase. More precisely by physical state we mean a point of the two dimensional 
sphere (called thc Bloch sphere) Sb '■= S 3 / ~ where the equivalence relation ~ is defined as follows: tp ~ ip' 
(where ip, ip' £ S 3 ) if and only if ip = exp (zí>)i/)', for some $ £ [0, 2ir[. The projection from S 3 to Sb is called 
Hopf projection and it is given explicitly in the next section. A particularly interesting transition is of course 
from the first to the second eigenstates of Ho (i.e., from \tpi\ 2 = 1 to \4>2\ 2 = !)■ 

Due to the presence of the drift, in this case the minimum time problem with bounded control and the 
minimum energy problem are different. In [20] the authors studied the minimum energy problem (in that case, 
optimal solutions can be expressed in terms of Elliptic functions), while here we minimize the time of transfer, 
with bounded field amplitudc: 

|fi(í)| < M, for every t £ [0,T], (8) 

where T is the time of the transition and M > represents the maximum amplitude available. This problem 
requircs complctely different techniques with respect to those used in [20]. 

Thanks to the reduction to a two dimensional problem (on the Bloch sphere) , this problem can be attacked 
with the techniques of optimal syntheses on 2-D manifolds developed by Sussmann, Bressan, Piccoli and the 
first author, see for instance [10, 17, 29, 35] and recently rewritten in [12]. We make a brief recali of these 
techniques in Appcndix A. 

1.3 The Control problem on the Bloch Sphere Sb 

An explicit Hopf projection from S 3 to Sb is given by: 

-2Re(V>í>2) 

n:( 02 \£S 3 dC 2 ^y=\ y 2 | = | 2Im(^0 2 ) | eS B cl 3 . (9) 
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State 2 



Figuro 1: The Bloch Spherc 



Notice that II maps the first eigenstate of Hq (i.c. iV'il 2 = 1) to the north pole P/v := (0, 0, 1) T of S^, and the 
second eigenstate (i.e. |t/í 2 | 2 = 1) to the south pole P s := (0,0, —1) T . 

After setting u(t) = fl(i)/M, the Schròdinger equation (6), (7) projeets to the following single input affinc 
system (clarified below, after normalizations) , 

V = F s {y)+uG s (y), \u\<l, where: (10) 

3 

yeS B :={( yi ,y 2 ,y 3 )£R\ = 1} 

3=1 

F s (y) :=kcos(a) | y x |, G s (y) := fcsin(a) | -y 3 |, (12) 
M 





a := aretan í — J 6 ]0, tt/2[, k := 2E/ cos(a) = 2yJ M 2 + E 2 . (13) 

Rcnmrk 2 (normalizations) In the following, to simplify the notations, we normalizc k — 1. This normal- 
ization corresponds to a reparametrization of the time. More preciscly, if T is the minimum timc to steer the 
state y to the statc y for the system with k = 1, the corresponding minimum time for the original system is 
T j(2\jM 2 + E 2 ). Sometimes we need also the original system (6), (7) on «S 3 , with the normalization made in 
this remark, i.e. the system 

,í§> =*(,>«(<), w het c fiW-Í* a (-* a ^i). (») 

We come back to the original value of k only in Section 3.3, where we compare our results with those of other 
authors. 

We refer to Figure 1. The vector fields Fs(y) and Gs(y) (that play the role respectively of Hq and H{) describe 
rotations respectively around the axes y 3 and y\. Lct us definc the vector fields corresponding to constant 
control ±1, 

Xf(y):=F s (y)±G s (y). (15) 

The parameter a €]0, 7r/2[ (that is the only parameter of the problem) is the angle between the axes of rotations 
of Fs and X§. The case a > 7r/4 (resp. a < 7f/4) corresponds to M > E (resp. M < E). 
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Definition 1 An admissible control u(.) for the system (10)-(13) is a measurable function u{.) : [a,b] — > [—1,1], 
while an admissible trajectory is a Lipschitz functions y(.) : [a,b] — > Sb satisfying (10) a.e. for some admissible 
control u(.). If y{.) is an admissible trajectory and u{.) the corresponding control, we say that (y(.),u(.)) is an 
admissible pair. 

For every y € Sb, our minimization problem is then to find the admissible pair steering the north pole to y in 
minimum time. More precisely 

Problem (P) Consider the control system (10)- (13). For every y € Sb, find an admissible pair (y(.), u(.)) 
defined on [0, T] such that y(0) = Pn, y{T) = y and y(.) is time optimal. 

In Optimal Control the problem (P) is known as the problem of computing the time optimal synthcsis for 
the system (10)-(13). For more elaborated definitions of optimal synthesis see Appendix A, or [12. 30] and 
references therein. 

Definition 2 (bang, singular for the problem (10)-(13)) A control u(.) : [a,b] — » [—1,1] is said to be 
a bang control if u(t) = +1 a.e. in [a,b] or u(t) = —1 a.e. in [a,b]. A control u(.) : [a,b] — > [—1,1] is said 
to be a singular control if u(t) = 0, a.e. in [a,b]. A finite concatenation of bang controls is called a bang-bang 
control. A switching time of u(.) is a time t € [et, b] such that, for every e > 0, u is not bang or singular on 
(í — e, t + e) fi [a, b] . A trajectory of the control system ( 27) is said a bang trajectory ( or arc ), singular trajectory 
(or arc), bang-bang trajectory, if it corresponds respectively to a bang control, singular control, bang-bang control. 
If t is a switching time, the corresponding point on the trajectory y(t) is called a switching point. 

Rcnmrk 3 The definitions of singular trajectory and control, given above are very specific to our problem 
(10)-(13). For the definition of singular trajectòries for more general systems sec Definition 8, Appendix A.l. 

In [11] it was proved that, for the same problem (10)-(13), but in which y € M.P 2 , for every couplc of points 
there exists a time optimal trajectory joining them. Moreover it was proved that every time optimal trajectory 
is a finite concatenation of bang and singular trajectòries. Repeating exactly the same arguments and recalling 
that S 2 is a double covering of M.P 2 , one easily gets the same result on Sg. More precisely we have: 

Proposition 1 For the problem (10)-(13), for each pair of points p and q belonging to Sb, there exists a time 
optimal trajectory joining p to q. Moreover every time optimal trajectory for the problem (10)-(13) is a finite 
concatenation of bang and singular trajectòries. 

Thus, the Fuller phenomenon (i.e. existenec of an optimal trajectory joining two points with an infmitc 
number of switchings in finite time) never oceurs. Notice that the previous proposition does not apply if a = 
or a = tt/2, since in these cases the controllability property is lost. 

1.4 Purpose of the paper 

Our aim is to study problem (P) for every possible value of the parameter a, giving a particular relief to the 
case in which y = Ps (i-e. to the optimal trajectory steering the north to the south pole). 

We will not be able to give a complete solution to the problem (P), without the help of numcrical simulations. 
However, thanks to the theory developed in [12] we give a satisfactory description of the optimal synthcsis. In 
the following we describe the main results and the structure of the paper. 

For a < 7r/4, every time optimal trajectory is bang-bang and in particular the corresponding control is 
periòdic, in the sense that for every fixed optimal trajectory the time between two consecutive switchings is 
constant. Moreover it tends to n as a goes to 0. For the original non normalized problem this means that for 
M/E « 1, the optimal control oscillates with frequeney of the order of the resonance frequeney cúr = 2E. 
In this case it is possible to give a satisfactory description of the optimal synthesis excluding a neighborhood 
of the south pole, in which we are able to computo the optimal synthcsis only numerically (such results were 
already present in [11] as we see below). 

On the other side. if a > 7r/4 the computation of the optimal synthcsis is simplcr since the number of 
switchings needed to cover the whole sphere is small (less or cqual than 2). In this case, for a big cnough. 
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we are also able to give the exact value of the time needed to cover the wholc sphere. However, therc is a 
new difficulty, namcly the presence of singular arcs. Moreover the qualitative shape of the synthesis is rathcr 
different if a is close to ir/4 or to ir/2. A relevant fact is that this synthesis contains a singularity (the so called 
(S,K) 3 ) that is predicted by the general theory (see [12], pag. 61 and 82), and was never observed out from 
ad hoc examples. 

The problem of finding explicitly the optimal trajectòries from the north pole P^ to the south pole P$, 
can be easily solved in the case a > ir/4 as a consequence of the construction of the time optimal synthesis. 
(Corning back to the original non normalizcd problem we also prové that fixed E, for M — > oo the time of 
transfer from P/v to Ps tends to zero.) 

For a < 7r/4 the problem is more complicated. However, using the symmetries of the problem, we are able 
to restrict the set S of candidatc optimal trajectòries reaching the south pole, to a set containing at most 8 
trajectòries (half starting with control +1 and half starting with control —1, and switching exactly at the same 
times). Thcse trajectòries are determined in terms of a parameter (the first switching time) that can be easily 
computed numerically solving suitable equations. Once these trajectòries are identified one can check by hands 
which are the optimal ones. 

The analysis can be pushed much forward. We also prové that the cardinality of S depends on the so called 
normalizcd remainder 



where [ . ] denotes the integer part. In particular, for a small, we prové that if R is close to zero then 5 contains 
exactly 8 trajectòries (and in particular there are four optimal trajectòries), while if R is close to 1 then S 
contains only 4 trajectòries (two of them arc optimal). The precise description of these facts is contained in 
Proposition 6. As a consequence, the qualitative shape of the time optimal synthesis presents different patterns, 
that cyclically alternate, in the non controllability limit a — ► 0, giving a partial proof of a conjecture formulated 
in a previous paper ([11]), that was supportcd by numerical simulations, see Rcmark 11. This is probably the 
most interesting byproduct of this paper. 

Finally we compare these results with some known results of Khaneja, Brockett and Glaser and with those 
obtained by controlling the magnètic field both on the x and y directions. 

The structure of the paper is as follows. In Section 2 we briefly resume the results of paper [11] which are 
connected to our problem and the conjectures formulated therein. The main results of the paper are described 
in Section 3, while the proofs are postponed to Appendix B. In Appendix A we recali the main tools of the 
theory of optimal synthesis. In Appendix C we determine the last point reached by trajectòries starting at Pjv 
and the time needed to cover the whole sphere. 

2 History of the problem and known facts 

The problem (P) (although with different purposes) was already partially studicd in [11], in the case a < 7f/4. 
In that paper the aim was to give an estimate on the maximum number of switching for time optimal trajectòries 
on 5*0(3) (problem first studied by Agrachev and Gamkrelidze in [2], using index theory). 

In [11] it has been proved that, for the problem (P) in the case a < 7r/4, every optimal trajectory is 
bang-bang. More preciscly, it was proved that in the case a < 7r/4, if y(.) is a time optimal trajectory starting 
at the north pole, then it should satisfy the following properties: 

i) y(-) is bang bang; 

ii) the duration s; of the first bang arc satisfies s; € [0,7r], 

iii) the time duration between two consecutive switchings is the same for all interior bang arcs (i.e. exeluding 

the first and the last bang) and it is the following function of Si defined in the interval [0,7r], 



One can immcdiatcly check that this function satisfies v(0) = v(ir) = ir and v(s;) > tt for every Sj s]0, 7r[, 




(17) 
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iiii) the time duration of the last arc is Sf <E [0, v(sí)}, 

Propcrties ï) — iiii) are illustrated in Figurc 2. Moreover, thanks to the analysis givcn in [11], one easily get 
(always in the case a < ir/4): 

v) the number of switchings N y of y(.) satisfies the following inequality 

+ 1 (18) 

Conditions i)-v) define a set of candidatc optimal trajectòries. The way in which these candidate optimal 
trajectòries cover the whole sphere is shown in the top of Figure 3. 

Consider the following curves, made by points where the control switches from +1 to —1 or viceversa, called 
switching curves, defined by induction 

Cf(s)=e x s v ^e x s es P N , Cl(s) = e x ^C^ 1 (s), (whexe e = ±1 and k = 2, N M - 1). (19) 
See the top of Figure 3. 

Even if the analysis made in [11] was sufficient to the purpose of giving a bound on the maximum number of 
switchings for time optimal trajectòries on 50(3), some qüestions remained unsolved. In particular qüestions 
about local optimality of the switching curves. Roughly speaking we say that a switching curve is locally 
optimal if it never "reflects" the trajectòries (see Figure 4 A). 1 When a family of trajectòries is reflected by a 
switching curve then local optimality is lost and some cut locus appear in the optimal synthcsis. 

Definition 3 A cut locus for the problem (P) is a set of points reached at the same time by two (or more) 
optimal trajectòries. A subset of a cut locus that is a connected C 1 manifold is called overlap curve. 

An example showing how a "reflection" on a switching curves generate a cut locus is portrayed in Figure 4 
B and C. More details are given later. In [11], the following qüestions remain unsolved: 



N y < N M :-- 



Question 1 Are the switching curves (7|, k = 1..., Nm — F locally optimal? More precisely, one would like to 
understand how the candidate optimal trajectòries described above are going to lose optimality. 

Question 2 What is the shape of the optimal synthcsis in a neighborhood of the south pole? 
Numerical simulations suggested some conjectures regarding the above qüestions. More precisely: 

Cl Define ki ast = [^^p] — 1- Then the curves C|(s), (k = 1, ...,N M — 1) are locally optimal if and only if 
k < hast- Notice that ki ast € {N M - 3, N M - 2}. 

1 More precisely consider a smooth switching curve C between two smooth vector field Y\ and Y2 on a smooth two dimensional 
manifold. Let C(s) be a smooth parametrization of C. We say that C is locally optimal if, for every s £ Dom(C), we have 
C(s) ^ aiYi(C(s)) + 02^2(^(5)), for every cei,ct2 s.t. «i«2 > 0. The points of a switching curve on which this relation is not 
satisfied are usually called "conjugatc points". See Figure 4. 
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Figurc 3: Synthcsis on thc sphcrc for a < ir/4 and conjectured shape in a ncighborhood of the south polc 
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Figure 4: Locally optimal switching curves and non locally optimal switching curves with the corresponding 
synthesis 
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Analyzing the evolution of the minimum time wave front in a ncighborhood of thc south-pole, it is reasonablc 
to conjecture that: 

C2 The shape of the optimal synthcsis in a ncighborhood of the south pole depends on the so called remaindcr 2 
r := 7r — 2a [j^] ■ Noticc that r belongs to thc interval [0, 2a[. More precisely, we conjecture that for 
a €]0, 7r/4[, there exist two positive numbers ol\ and ai such that < a± < a < ai < 2a and: 

CASE A: r e.]a 2 ,2a[. The switching curve t glues to an overlap curve that passes through the 

origin (Fig. 3, Case A). 

CASE B: re [a±, a2\- The switching curve C £ Nm _ 1 is not reached by optimal trajectòries in the interval 
JO, 7rJ. At thc point x (0) an overlap curve starts and passes through the origin. 

CASE C: r G]0, ctx[. The situation is more complicated and it is depicted in the bottom of Fig. 3, 
Casc C. 

For r = 0, the situation is the same as in CASE A, but for the switching curve starting at C^, M _ 2 (0). 



3 Main Results 

We give here a brief description of thc main results of the paper. The corresponding proofs are given in 
Appcndix B. From now on we use the following conventions. 

Remark 4 (notation) The letter B refers to a bang trajectory and the letter S refers to a singular trajectory. 
A concatenation of bang and singular trajectòries is labclcd by thc corresponding letter sequence, written in 
order from left to right. Sometimcs, we use a subscript to indicate the time duration of a trajectory so that we 
use B t to refer to a bang trajectory defined on an interval of length t and, similarly, St for a singular trajectory 
defined on an interval of length t. Moreover we indicate by 7+ (resp. 7~) the trajectory of (10)-(13) starting 
at the north pole at time zero and corresponding to control u = 1 (resp. u = — 1). Notice that 7 ± are defined 
for every time, and are periòdic. Finally we use the following subsets of S^: the circle of equation y 3 = called 
equator, the set 7/3 > 0, called north hcmisphere and the set 7/3 < 0, called south hcmisphere. 



3.1 Optimal synthesis for ot > 7r/4 

In this section we describe the time optimal synthesis for a > 7r/4. We divide Sb in 8 open regions called 
íi^, Çlf, ^nasty anc ^ m ^ arcs ( see Defmition 4, and Figure 5). For every point y £ Sb \ (^n as t y U ^nasty)> 
Theorem 1 gives the optimal trajectòries reaching y. 

Unlike the a < 7r/4 case, here it is possible to detect the presence of singular trajectòries that are optimal, 
and also of cut loci (even not only in a ncighborhood of the south pole). 

The region í^asty (and similarly íín asty ) is more dimcult to analyzc. It contains a cut locus that should be 
determined numcrically. Even if we are not able to provide an analytic characterization of this locus, we are 
able to prové the following. 

i) a = arcsin(l/ v^2) is a bifurcation point for the optimal synthesis i.e. the qualitative shape is different if 

a € [7r/4, arcsin(l/v / 2)[ (called Case 1) or a <G [arcsin(l/%/2), 7r/2[ (called Case 2). More precisely, from 
the point D + := 7 + (7r), in Case 1 it starts an optimal switching curve, whilc in Case 2 it starts an 
overlap curve (see Proposition 3). The situation in íí~ asty is symmetric. 

ii) The south pole belongs to the cut locus and it is reached exactly by four optimal trajectòries (see Proposi- 

tion 2). 

Numerical computations show that in Case 2, thc cut locus in í^n asty is an overlap curve connecting D + 
with the south pole, while in Case 1, the switching curve starting from D + loscs local optimality at a point 
of íín as ty and conneets to an overlap curve which reaches the south pole (see Figure 6). Remark 9 explains 

2 Notice that r = 2aR, where R has been defined in Formula (16). In conjecture C2, we use the remainder r, to keep the same 
notation of [11]. 
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that in Case 2 it is not necessary to compute the cut locus lying in í^asty ^° ^he expression of the optimal 
trajectory connecting Pjy to a point of í^nasty The situation in í^nasty ^ s symmetric. 

Let us start with the description of the optimal synthesis in Sb \ (^íasty ^ ^ïïasty)- Even if Definition 4 
and Theorem 1 look complicated, the shape of the optimal synthesis is quite simple as it is shown in Figure 6. 

Definition 4 According to Figure 5, let us define the following curves on Sb- 

• Let t\ be the first time at which 7 + interseets the equator and let A + := j + (ti) (notice that t\ = 
7r — arccos(cot 2 (a)) / ). Define PnA+ = Supp (7 + |[o.t 1 ])· 

• Let Ç~ be the trajectory corresponding to control —1, starting at time zero from A + . Let ti be the first 
positive time at which Ç~ interseets the equator (notice thatt 2 = 2 arccos(cot 2 (a)) ). Define B + := £"(£2) 



and AB+ = Supp (£ \[o,t 2 ])- 



• Let + = (1,0,0). Define AO + (resp. OB + ) as the support of the trajectory corresponding to control 
zero, starting at A + (resp. + ) and ending at + (resp. B + ). 



• Recali that D + = 7 + (7r), and define AD + = Supp (7~*"|[ti,7r])j DB + = Supp \f^\w,t3\) , where Í3 is the 
second intersection time of j + with the equator (notice that Í3 = ir + arccos(cot 2 (o!)) = t\ + ti)- 



• Let BPg the support of the trajectory corresponding to control — 1, starting at B + and ending at the 
south pole. 



• Let DPg the connected subset of the meridian yi = 0, lying in the south hemisphere and connecting the 
point D + to the south pole. 



Similarly define A~ , B~ , 0~ , D~ , P N A~, AB~ , AO~ , OB~ , AD~ , DB- , BP S , DP S . 

According to Figure 5 define Í7 X , . . . , fij, Qnasty as ^ e °P en connected components of the open set obtained 
subtracting from Sb all the arcs defined above. 

The following theorem holds for every a G]7r/4, 7r/2[. For the particular valuc a = 7r/4 the claims of the theorem 
must be modihed. Such changes are reported in Remark 5. 

Theorem 1 LefYy be the set of time optimal trajectòries steering the north pole to y. We have the following: 



Tl. // y G PnA + then 'Yy is made by a unique trajectory corresponding to control +1 of the form Bt, with 
t < t x . 
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T2. If y G AB + \ B + then "y y is made by a unique trajectory of the form B tl B t (wíth the first bang corre- 
sponding to control +1). 



T3. If y G AO + then "Yy is made by a unique trajectory of the form Bt 1 S s (with the first bang corresponding 
to control +1). 



T4. If y G 0B + \ + then "Yy is made by two trajectòries of the form B tl S s B t , both starting with control +1 
and ending respectively with control +1 and — 1. These two trajectòries have the same vàlues of s > 
and t > 0. 



T5. If y G AD + then y g ós made by a unique trajectory corresponding to control +1 of the form B t , with 
t G [íi,7r]. 



T6. If y £ DB + \ B + then / Yg is made by a unique trajectory corresponding to control +1 qf í/ie /orm Bt, 
wraí/i í G Í7T, Í3[. 



T7. I/y G -B-Pg then "^y is made by two trajectòries respectively of the form Bt 1 B t and B ta B t -t 2 an d starting 
with control +1. 



T8. Ifyeíl+U (DP S \ P s ), then fy is made by a unique trajectory of the form B t B t ' , with < t < t\ and 
the first bang corresponding to control +1. 

T9. If y G , then 'Yy is made by a unique trajectory of the form B tl S s B t , with s > 0, the first bang arc 
and the last bang arc corresponding respectively to control +1 and — 1. 

T10. If y G Qïj", then *y y is made by a unique trajectory of the form B tl S s B t , with s > and both bang arcs 
corresponding to control +1. 

Tll. If y = Ps then 'V y is made by the four trajectòries of the form B tl B t3 and B ta B tl . 

T12. If y G Qnasty ^ en every trajectory of^g is bang-bang with at most two switchings. 

If y belongs to one of the remaining sets defined above, the description of the optimal strategy is analogous, by 
symmetry. 

Rcmark 5 In the case a = 7r/4 some changes in the previous statement are required. In particular the points 
A + ,B + ,0 + and D + coincide (also the points A~ ,B~ ,0~ and D~ coincide) and, consequently, there are no 
optimal trajectòries containing singular arcs. Another immediate consequence of this fact is that there are only 
two optimal trajectòries reaching the south pole, of the form B^B^. 



Rcmark 6 Notice that every point of OB + \ + , OB \ O , BPg , BP S is reached by more than one optimal 
trajectory, i.c. it belongs to the cut locus. Other points of the cut locus can be identified numerically in íln asty 
and í^nasty as explained in the next section. 

Rcmark 7 In Theorem 1 we do not specify all the durations of the bang arcs. However the missing ones can 
be obtained simply by following the switching strategy backwards. 

Remark 8 Note that the region reached by optimal trajectòries containing a singular arc nfunf uao^uob ± 
become bigger and bigger as a tends to n/2. Moreover, in this limit, since the modulus of the drift Fs becomes 
smaller and smaller, the time needed to cover such region tends to infinity. Notice however that the time needed 
to reach Ps is always 2tt. The time needed to reach every point of the sphere for a big enough, and the last 
point reached by an optimal trajectory containing a singular arc, can be computed explicitly. This is done in 
Appendix C. 

Since the case y = Ps is important also for the determination of the cut locus in íín asty U ^n as t yI h is 
reported in the next section as a separate proposition (see Proposition 2). 
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Case 1 Case 2 



Figuro 6: Optimal synthesis for a = 7r/3 and a slightly larger than 7r/4. 



3.1.1 The time optimal synthesis in í^n asty and optimal trajectòries reaching Ps for a > ir/4 

From next proposition. Tll of Theorem 1 follows. Morc prcciscly Proposition 2 shows that in the case a > 7r/4, 
thcre are exactly four optimal trajectòries steering P/v to Ps, and it characterizes them. As a consequence, the 
south pole belongs to the cut locus. 



Proposition 2 Consider the control system (10)-(13), and assume a > 7r/4. Then the optimal trajectòries 
steering the north pole to the south pole are bang-bang with only one switching. More precisely they are the four 
trajectòries corresponding to the four controls 

( i)_ / u = i, te[o,h] (2) _ r i, íg[o,í 3 ] „(3)_j -i, *G[o,ti] (4) _ r -i, íg[o,í 3 ] 
-\ u =-i, t€}t U T], u \ -i, te]f 3 ,n \i, t€]íi,n íg]í 3 ,t] 

where t\ and Í3 are defined in Definition 4, cind T = 2ir. 



One can easily check that the switchings described in Proposition 2 oceur on the equator (y^ = 0). 

The following proposition describes the optimal synthesis in f^asty m a neighborhood of the points D ± 
and the bifurcation oceurring at a = arcsin(l/v2)- 

Proposition 3 Let a > 7r/4. /na neighborhood of the point D + in Qnasty there exists a switching curve 
starting at D + of the form é"^ x s e sX s Pm ■ If a> 7r/4 thís curve is tangent to the equator at D + . Moreover if 
a < arcsin(l/ \/2) ( above called Case 1 ) then the switching curve is optimal near D + , while if a > arcsin(l/v^2) 
( above called Case 2 ) then the switching curve is not locally optimal near D + and an overlap curve starts at 
the point D + . A symmetric result holds in a neighborhood of D~ in Qnasty 

The region í^nasty contains a cut locus that should be determined numerically. In Case 2, numerical 
simulations show that the switching curve starting at D + is never optimal, i.e. every point of í^asty ^ s 
reached by an optimal trajectory of the form e +e Pjv, with s G]0, íi[ or an optimal trajectory of the form 
e tx -e sX +P N , with s g]7t,í 3 [. 
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Figurc 7: Graph of v(.) when a = n/6 (figure A). Graph of the functions T and Q when a = 0.13 (figure B 
and C) 



Rcmark 9 Notice however that, in Case 2, givcn a point y G ^nasty to find thc timc optimal trajectory 
reaching y, it is not necessary to compute the cut locus. Indeed it is sufficient to compare the final times, 
corresponding to the two switching strategies given above, and to chose the quickest one. The situation in 
^nasty is symmetric. 

In Case 1, the situation is morè complicated. The switching curve dcscribed by Proposition 3 has the 
expression C+(s) = e x s< s )e x s s P N , s S]0,íi[ where the function v(.) is given by the same formula of the 
q < n/A case, i.e. v(s) = n + 2 arctan [(sins)/ (cos s + cot 2 a)] . (To verify such formula it is enough to repeat 
the computations done in [11].) As described by Proposition 3, this switching curve is optimal near D + and 
numerical simulations show that there exists s G]0, íi[ such that there is an optimal trajectory switching on 
C^"(s) if and only if s G [0, s[, and an overlap curve connecting C^~(s) to the south polc appears. The optimal 
synthcsis for Case 1 and Case 2 is depicted in Figure 6. 

3.2 Optimal trajectòries reaching the south pole for a < n/A 

In this section we characterize the time optimal trajectòries reaching the south pole, in the case a < ■k/A. 
This characterization is more complicated with respect to the case a > n/A, due to the fact that the optimal 
trajectòries have many switchings. The time optimal synthesis for a < n/A was already (partially) studied in 
[11] and it has been described in Section 2. 

From conditions i) — iiii) in Section 2, we know that every optimal trajectory starting at the north pole has 
the form B Si B v r 3i \ ■ ■ ■ B v r 3i \B s , where the function v(sí) is given by formula (17). (In the following we do not 

specify if the first bang corresponds to control +1 or —1, since, as a consequence of the symmetries of thc 
problcm, if u(t) is an optimal control steering the north pole to the south pole, —u(t) steers the north pole to 
the south pole as well.) It remains to identify one or more vàlues of Sj, Sf and the corresponding number of 
switchings n for this trajectory to reach the south pole. 

Notice that t = arccos(— tan 2 (a)) is the maximum of the function i?(.) on the interval [0, n], v(.) is increasing 
on [0,t\ and decreasing on [t,n] and v(0) = v(n) = n. Thcn, given s G [0,n] such that s^í, there is a unique 
solution s*(s) G [0, 7r], s*(s) ^ s, to the equation v(s*) = v(s). Thc function s*(.) is extended to the whole 
interval [0, n] setting s*(t) = ï (see Figurc 7 A). Thanks to the symmetries of the problcm, we prové that 
if a < n/A, Sf is equal either to s, or to s*(sí). This fact is described by Lemma 4 stated and proved in 
Appcndix B. 

The following two propositions describe how to identify candidatc triples (sj, s/, n) for which the correspond- 
ing trajectory steers the north pole to the south pole in minimum time. We say that a bang-bang trajectory, 
solution of the system (10)-(13), is a candidate optimal trajectory if it is an extrcmal trajectory for problcm 
(P) reaching the south pole and it has a number n of switchings satisfying n < Nm (defined in Formula (18)). 
From Lemma 4, there are two kinds of candidate optimal trajectòries: 

• Sf = s*(sí), callcd TYPE-l-candidatc optimal trajectòries 

• Sf = Si called TYPE-2-candidate optimal trajectòries 
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Dcfinc the following functions, whose geomètric meaning is clarificd in Appendix B.2: 

6(s) = 2 arecos íshi 2 (^) cos(2a) - cos 2 J ) (20) 



[3(s) = 2 arccos(sin(a) cos(a)(l — cos(s))) (21) 
Proposition 4 (TYPE-l-trajectories) Fixcd a < 7r/4, the equation for the couple (s,n) £ [0,ir] x N: 

^):^=", (22) 

Las either two or zero solutions. More precisely if (s,n) is a solution to equation (22), then (s*(s),n) is the 

second one. The trajectòries B s B v r s \ ■ ■ ■ B v r s \ -B s »( s ) and B g »í g \ B v r s -\ ■ ■ ■ B v r s \ B s are the TYPE-l-candidate 

- . ' . ' 

n—t n—1 

opthnal trajectòries. 

Proposition 5 (TYPE-2-trajectories) Fixcd ot < 7r/4, the equation for the couple (s,n) £ [0, ir] x N: 

S(S): =W +1 =„, (23 ) 

has exactly two solutions. More precisely these solutions have the form (si,n), (s2,n + 1). The trajectòries 
B Sl B v ( Sl ) ■ ■ ■ B v r Sl \ B Sl and B S2 B v t S2 \ ■ ■ ■ B„( S2 ) B S2 are the TYPE-2-candidate optimal trajectòries. 

S v ' V v ' 

n—1 n 

In Figure 7 B and C the graphs of the functions (22) and (23) are drawn for a particular value of a, namely 
et = 0.13. Propositions 4 and 5 select a set of (possibly coinciding) 4 or 8 candidate optimal trajectòries (half 
of them starting with control +1 and the other half with control —1) corresponding to triples (sj, Sf,n). Such 
triples can be easily computed numcrically solving cquations (22) and (23). Then the optimal trajectòries can 
be selected by comparing the times needed to reach the south pole for each of the candidate optimal trajectory. 
Notice that there are at least two optimal trajectòries steering the north to the south pole (one starting with 
control +1 and the other with control —1). 

If 7r/(2a) is an integer number n, then TYPE-1 candidate optimal trajectòries coincide with the TYPE-2 
candidate optimal trajectòries of the form B^ B^.·.B^ B^. The remaining trajectòries of TYPE-2 are of the 

n— 2 

form B s B v ( s y..B v ( s ) B s for some s e]0, 7r[. Otherwise if tt/ (2a) is not an integer number, define: 



7T . 7T 

í — 1, and the normalizcd remainder R := — 
L 2a 2a 



7T 

2a 



G [0,1[ 



where [.] denotes the integer part. The following proposition determines precisely the time optimal trajectòries 
for particular vàlues of the parameter R: 

Proposition 6 For m large enough there exist ri(m) < r^ijn) g]0, 1[ such that: 

A. if R £]0,ri(m)] then equation (22) admits exactly two solutions that are both optimal, while TYPE-2 can- 

didate optimal trajectòries are not. 

B. z/R £\r\(m),r2(m)[, then equation (22) admits two solutions, that are not optimal. 

C. i/R £\r2(m), 1[ then equation (22) does not admit any solution. Moreover ^(rn) — > for m — > oo. 

Rcnmrk 10 The function r2(m) can be determined explicitly (see Appendix B.2.1), while for r\(m) we are just 
able to prové the existence, and we conjecture that it can be taken equal to r2(m). 
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Rcnmrk 11 An important conscqucnce of Proposition 6 is that for a small, the number of optimal trajectòries 
reaching the south pole is not fixed with respect to a. Indeed such number altcrnates as a — > 0, according to 
Proposition 6: in particular it is equal to 4 if R <E]0,ri(m)] and it is equal to 2 if R e]r 2 (m), 1[U{0}. This is 
enough to conclude that also the qualitative shape of the optimal synthesis in a neighborhood of the south pole 
alternates giving a partial proof to the conjecture C2 of Section 2 (originally stated in [11]). In particular it is 
a proof of the first assertion (on the dependence of the synthesis on the remainder r = 2aR). Moreover notice 
that the results of Proposition 6 perfectly fit with all the other statements of conjecture C2 with r2(m) playing 
the role of ai/(2a). One can apply the definition of locally equivalent syntheses given in [12] (see Dcfinition 
32, pag. 59), to make rigorous the statement that the qualitative shape of the optimal synthesis changes with 
a. 

Using the previous analysis one can easily show the following result (of which we skip the proof): 

Proposition 7 If N is the number of swítchings of an optimal trajectory joining the north to the south pole, 
thcn 

7T 7T 

1< N < l·l. 

2a ~ 2a 

Using these inequalities and the fact that, for a < n/6, the function 2s + ( 1 ) v(s) is inercasing on 



[0, 7r], one can give a rough estimate of the time needed to reach the south pole: 



2o 



Proposition 8 The total time T of an optimal trajectory joining the north to the south pole satishes the 
inequalities: 

TT 2 TT 2 

2n <T < hTT. 

2a 2a 

3.3 Comparison with results in the rotating wave approximation and with [27] 



In this section we come back to the original value of k i.e. k = 2E/ cos(a) = 2V M 2 + E 2 , and we compare the 
time necessary to steer the state one to the state two for our model and the model (4), described in Remark 1, 
in which we control the magnètic field both along the x and y direction, or we consider a two-level molecule 
in the rotating wave approximation. We recali that — E, E are the energy levels and M is the bound on the 
control. For our model, the time of transfer T satisfies: 



• for ol > tt/4 (i.e. for M > E) thcn T = 2n/k = n/yjM 2 + E 2 ; 

• for a < 7r/4 (i.e. for M < E) thcn T is estimated by \ ( ^ 2n) < T < i + n 

k \2a J k y 2a 

On the other hand, for the model (4), the time of transfer is Te = n/(2M) (cfr. Remark 1). Fixed E = 1, 
in Figure 8 A the times T and Te as function of M are compared. Notice that although Te is bigger than 
the lower estimate of T in some interval, we always have Te < T. This is due to the fact that the admissible 
velocities of our model are a subset of the admissible velocities of the model (4). 

Notice that, fixed E = 1, for M -> we have T - tt 2 /(4M) = (n/2)T c , while for M -> oo, we have 
T - n/M = 2T C . 

Remark 12 For M << E (i.e. for a small) the difference between two switching times is v(s)/k ~ n/(2E). 
It follows that a time optimal trajectory connecting the north to the south pole (in the interval between the 
first and the last bang) is periòdic with period P ~ tt/E i.e. with a frequeney of the order of the resonance 
frequeney lur = 2E (see Figure 8 B). On the other side if M > E thcn the time optimal trajectory connecting 
the north with the south pole is the concatenation of two pulscs. Notice that if M » E, the time of transfer 
is of the order of tt/M and therefore tends to zero as M — * oo. It is interesting to compare this result with 
a result of Khaneja, Brockett and Glaser, for a two level system, but with no bound on controls (see [27]). 
They estimate the infimum time to reach every point of the whole group SU(2) in n/E. On the other side, in 
Appcndix C it is proved that the time needed to cover the whole sphcrc S B = SU(2)/S 1 goes to n/(4E) as 



16 



Re(Q) 




Figurc 8: A. Estimatc of thc minimum time to reach the state two and comparison with the time needed with 
two controls or in thc rotating wave approximation B. Comparison between the optimal strategy for our 
system and in the rotating wave approximation 



M goes to infinity (however this does not contradict the fact that the state two can be reached in an arbitrary 
small time, as we discussed above). 

Notice that our optimal control has the same form of the control computed in [27] i.e. a pulse (bang) 
followed by an evolution with the drift (singular) followed by a pulse (bang). 

A An overview on Optimal Synthesis on 2-D Manifolds 

In this scction wc briefly recali the theory of optimal synthcscs on 2-D manifolds for system of the kind 
y = F(y) +uG(y), \u\ < 1, developed by Sussmann, Bressan, Piccoli and the first author in [10, 17, 29, 35] and 
recently rewrittcn in [12]. This appendix is written to be as much self-consistent as possible. 

For every coordinate chart on the manifold it is possible to introduce the following three functions: 

A A (y) := Det(F(y),G(y))=F 1 (y)G 2 (y)-F 2 (y)G 1 (y), (24) 
A B (y) := Det(G(y),[F,G](y))=G 1 (y)[F,G} 2 (y)-G 2 (y)[F,G] 1 (y), (25) 
fs(y) := -A B {y)/A A (y). (26) 

The sets A^ 1 (0), A^ 1 (0) of zeroes of A^, A B are respectively the set of points where F and G arc paral·lel, and 
the set of points wherc G is paral·lel to [F, G] . Thcsc loci arc fundamcntal in thc construction of thc optimal 
synthesis. In fact, assuming that they are smooth embedded one dimensional submanifold of M wc have the 
following: 

• in cach connected region of M \ (A^ 1 (0) U A^ 1 (0)) , every extremal trajectory is bang-bang with at most 
one switching. Moreover, for every switching of the extremal trajectory the value of the control passes 
from -1 to +1 if fs > and from +1 to -1 if f s < 0; 

• thc support of singular trajectòries (that are trajectòries for which thc switching function idcntically 
vanishes, see Defmition 7 below) is always contained in the set A^ 1 (0); 

• a trajectory not switching on the set of zeroes of G is an abnormal extremal (i.e. a trajectory for which 
the Hamiltonian given by the PMP vanishing, see below) if and only if it switches on the locus A^ 1 (0). 

Then the synthesis is built recursively on the number of switchings of extremal trajectòries, canceling at 
each step the non optimal trajectòries (see [12], Chapter 1). 

Rcmark 13 Notice that, although the functions Aa and Ab depend on the coordinate chart, the sets A^ 1 (0), 
A^ 1 (0) and the function fs do not, i.e. they are intrinsic objeets of the control equation y = F(y) + uG(y). 
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A.l Basic Definitions and PMP on a n-dimensional Manifold 

In this section we define our optimization problem, we state the Pontryagin Maximum Principle, and we givc 
some bàsic definitions in the more general case of a n-dimcnsional manifold. We do this, since in Appcndix B.l 
we stated some result for the original problem (14), on S 3 ~ SU(2). 

Problem (Q) Consider the control system: 

y = F(y)+uG(y), y G M, \u\ < 1, (27) 

where: 

(HO) M is a smooth n-dimensional manifold. The vector fields F(y) and G{y) are C°°. 

We are interested in the problem of reaching every point of M in minimum time from a point yo G M. 

Definition 5 An admissible control u(.) for the system (27) is a measurable function u(.) : [a, b] — > [—1,1], 
while an admissible trajectory is a Lipschitz functions y{.) : [a, ()]-»¥ satisfying y(t) = F(y(t)) + u(t)G(y(t)) 
a.e. for some admissible control u(.) 

In the following we assume that the control system is complctc i.c. for every measurable control function 
u(.) : [a,b] —> [—1,1] and every initial state y, there exists a trajectory y(.) corresponding to «(.), which is 
defined on the whole interval [a, b] and satisfies y(a) = y. For us a solution to the problem (Q) is an optimal 
synthesis. For a more elaborated definition of synthesis, see Section A.3. 

Definition 6 (Optimal Synthesis) An optimal synthesis for the problem (Q) is a collection of time optimal 
trajectòries T = {y Vl (.) : [0,b yi ] h-> M, yi € M : y Vl (Q) = y , y Vl (b vl ) = yi}. 

The kcy tool is the PMP (see [3, 12, 25]). 

Theorem (Pontryagin Maximum Principle for the problem (Q)) Consider the control system (27) 
subject to (HO). Define for every (y, \,u) € T*M x [—1, 1] the function 

Jf(y,X,u) := < X,F(y) > +u < X,G(y) > . 

If the couple (y(.),u(.)) : [0, T] — > AI x [—1,1] is time optimal then there exist a never vanishing Lipschitz 
continuous covector A(.) : t £ [0, T] i— > A(í) G T*^M and a constant Ao < such that for a.e. t G [0,T]: 

m = -g^(y(t),Kt),u(t)), 

r) W 

ii) A(í) = — _(y( t ),A(í),u(t)) = -A(t)(VF(y(t)) + «(í)VG(ï(í))), 

iii) Jf(y(t),X(t), u(t)) = Jf M (y(t), X(t)) where Jíf M (v, X) := max{Jf (y, X,u) : w G [-1, 1]}, 
iiii) Mk(y(t), A(í)) + A = 0. 

Remark 14 The PMP is just a necessary condition for optimality. A trajectory y(.) (resp. a couple (y(-), A(.))) 
satisfying the conditions given by the PMP is said to be an extremal (resp. an extremal pair). An extremal 
corresponding to Ao = is said to be an abnormal extremal , otherwise we call it a normal extremal . 

We are now interested in determining the extremal trajectòries satisfying the conditions given by the PMP. A 
key role is played by the following: 

Definition 7 (switching function) Let (?/(.), A(.)) be an extremal pair. The corresponding switching function 
is defined as <$>(£) :=< A(í),G(y(í)) >. 

Notice that <j)(.) is continuously diffcrcntiable (indeed 4>(t) =< A(í), [F, G](y(t)) >, that is continuous). 
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Definition 8 (bang, singular) Letj, defined in [a,b], be an extremat trajectory andu(.) : [a,b] — > [—1,1] the 
corresponding control. We say that u(.) is a bang control if u(t) = +1 a.e. in [a,b] or u(t) = —1 a.e. in [a,b]. 
We say that u(.) is singular if the corresponding switching function <f>(t) = in [a,b]. A finite concatenation 
of bang controls is called a bang-bang control. A switching time of u(.) is a time t G [a,b] such that, for every 
e > 0, u is not bang or singular on (t — e, t + e) n [a, b] . An extremat trajectory of the control system ( 27) is said 
a bang extremat, singular extremat, bang-bang extremat respectively, if it corresponds to a bang control, singular 
control, bang-bang control respectively. If t is a switching time, the corresponding point on the trajectory y(t) 
is called a switching point. 

The switching function is important because it determines where the controls may switch. In fact, using the 
PMP, one easily gets: 

Proposition 9 A necessary condition for a time t to be a switching is that 4>(t) = 0. Therefore, on any interval 
where <fi has no zeroes ( respectively finitely many zeroes ), the corresponding control is bang ( respectively bang- 
bang). In particular, <fi > (resp <fi < 0) on [a,b] implies u = 1 (resp. u = — 1) a.e. on [a,b]. On the other 
hand, if (j> has a zero at t and 4>(t) is different from zero, then t is an isolated switching. 



A.2 More on singular extremals and predicting switchings for 2-D systems 

Now we come back to the case in which M is two dimcnsional. In this Scction we compute the control 
corresponding to singular extremals and we would likc to prcdict which kind of switchings can happcn, using 
properties of the vector fields F and G. The following two lcmmas illustrate the role of the functions A^ 1 (0), 
A^ 1 (0) in rclation with singular and abnormal extremals. The proofs can be found in [10, 12, 29]. 

Lemma 1 Let y(.) be an extremat trajectory that is singular in [a,b] C Dom{y{.)). Then y(.)|r a ,6] corresponds 
to the so called singular control <p(y(t)), where: 

VA B (y)-F(y) 

^ = -vmW)' (28) 

with Aa and Ab defined in Eqs. (24) and (25). Moreover, on Supp(y(.)), ip(y) is always well-defined and its 
absolute value is less than or equal to one. Finally Supp(y(.)\\ a u) C A^ 1 (0). 

Lemma 2 Let y(.) be a bang-bang extremat for the control problem (27), to G Dom(y{.)) be a time such that 
4>{to) = and G(y(to)) ^ 0. Then, the following conditions are equivalent: i) y(.) is an abnormal extremat; ii) 
y(t ) G A^ 1 (0); iii) y(t) G A^ 1 (0) ; for every time t G Dom(y(.)) such that <f>(t) = 0. 

The following lemma describes what happens when A^ and A# are different from zero. 

Lemma 3 Let íl C M be an open set such that il PI (A^ 1 (0) U A^ 3 1 (0)) = 0. Then all connected components 
of Supp(y(.)) n Í2, where y(.) is an extremat trajectory of (27), are bang-bang with at most one switching. 
Moreover, if fs > throughout Í2, then y(-)\n is associated to a constant control equal to +1 or —1 or has a 
switching from —1 to +1. If fs < throughout Í2, then y{-)\çi is associated to a constant control equal to +1 
or —1 or has a switching from +1 to — 1. 



A.3 Prame Curves and Frame Points 

For the problem (Q), under genèric conditions on the vector fields F and G, one can make the completo 
classification of synthcsis singularitics, stable synthcsis, singularities of the minimum time wave fronts. 

In the following, for sake of completeness, we recali the main results on existence of an optimal synthcsis 
and on classification of synthesis singularitics obtained in [28, 29] (sec also [12]). In [28], it was proved that the 
control system (27), under genèric conditions on F and G (with the additional assumption F(yo) = 0) admits 
a time optimal regular synthesis in finite time T, starting from yo. By genèric conditions, we mean conditions 
verified on an open and dense subset of the set of C°° vector fields endowed with the C 3 topology (see [12], 



19 



FCs ofkindY 
(X similar) 



FCofkindC FCofkindS FCofkindK 



Figure 9: Frame Curves of the optimal synthesis, under genèric conditions 



formula 2.6 pp. 39). To define what we mean by regular synthesis, we first nccd to introduce the concept of 
reachable set and of stratification of the reachable set. We call reachable set in time T > 0, the set: 

K(T) := {y 6 M : 3 b y G [0,T] and a trajcctory 

-f y : [0, b y ] -> M of (27) such that 7 y (0) = y 0l J y (b y ) = y}. 

Thcn we necd the dcfmition of stratification of 1Z(T), that roughly speaking is a partition in manifolds of 
diffcrcnt dimensions. 

Definition 9 (stratification) A stratification oflZ(T), T > 0, is a finite collection {Mi} of connected embed- 
dedC 1 submanifolds ofM, called strata, such thatlZ(T) = UkMk, and the following holds. If MjDClos(Mk) ^ 
with j 7^ k thcn Mj C Clos(Mk) and dim(Mj) < dim(Mk). 

Then a time optimal regular synthesis is defined by: i) a family of time optimal trajectòries T ~ {j y : [0, b y ] — > 
M, y G 1Z(T) : J y (0) = yo, ly(b y ) = y} such that if -f y G T and y = j y M for some t £ [0,b y ], then 
"/y = 7j/|[o,í]; ii) a stratification of 1Z(T) such that the optimal trajectòries of T can be obtained from a fcedback 
u{y) satisfying: 

• on strata of dimension 2, u(y) = ±1, 

• on strata of dimension 1, called frame curves (FC for short), u(y) = ±1 or u{y) = <p(y), where f{y) is 
defined by (28). 

The strata of dimension are called frame points (FP). Every FP is an intersection of two FCs. In [29] (sec 
also [12]), it is provided a complete classification of all types of FPs and FCs, under genèric conditions. All the 
possible FCs are: 

• FCs of kind Y (resp. X), corresponding to subsets of the trajectòries j + (resp. 7") defined as the 
trajectory starting at yo with constant control +1 (resp. constant control —1); 

• FCs of kind C, called switching curves, i.e. curves made of switching points; 

• FCs of kind S, i.e. singular extrcmals; 

• FCs of kind K , called overlaps and reached optimally by two trajectòries coming from different directions; 

• FCs which are arcs of optimal trajectòries starting at FPs. These trajectòries "transport" special infor- 
mation. 

The FCs of kind Y, C, S, K are depicted in Fig. 9. There are eighteen topological equivalence classes of FPs. 
A detailed description can be found in [9, 12, 29]. 
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B Proof of the Main Results 



In this section we give the proof of our main results. We start with a lemma, stating a property of optimal 
trajectòries, that is a consequence of the symmetries of the problem. It is used to identify the time optimal 
trajectòries steering the north to the south pole both for a > 7r/4 and a < ir/4. 

Lemma 4 Let a G]0, 7r/2[. Every optimal bang-bang trajectory, connecting the north to the south pole, with 
more than one switching is such that v(sí) = v(sf) where Si is the first switching time, Sf is the time needed 
to steer the last switching point to the south pole and u(sj) is the time between two consecutive switchings. 

Proof of Lemma 4. Considcr the problem of connecting Ps with P/v in minimum time for the system 
z = F' s (z) + uG' s (z) where z g S 2 and F' s (z) = —Fs(z), G' s (z) — —Gs(z). The trajectòries of this system 
coincide with those of the system (10)-(13), but the velocity is reversed. Thcrefore the optimal trajectòries for 
the new problem coincide with the optimal ones for the system (10)-(13) connecting P/v to Ps, and the time 
between two switchings is the same. Since performing the change of coordinates (z±, 22,^3) — * (2/1,2/2,2/3) = 
(— Z\, Z2, — Z3), the new problem becomes exactly the original problem, we deduce that, if we have more than 
one switching, it must be v(sí) = v(sf). I 



B.l Time Optimal Synthesis for the two Level Quantum System for a > 7r/4 

In this section, we apply the theory of optimal syntheses on 2-D manifolds recalled in Appendix A, to the 
system (10)-(13). Our aim is to describe the time optimal synthesis for a > 7r/4, i.c. to prové Thcorcm 1 
and Propositions 2 and 3. First we state some general results, holding for a g]0,7t/2[, regarding time optimal 
trajectòries of the system (14), on S 3 ~ SU(2), analogous to those obtained in [11] for SO(3) (in particular 
the proofs can be repeated using the same arguments). 

B.l.l General results on S 3 

In this section a s]0, 7r/2[. The first proposition states that singular extrcmals, defined as extremals for which 
the switching function vanishcs (sec Dcfinitions 7 and 8) correspond to zero control. This fact is very specific 
for our problem. 

Proposition 10 For the normalized minimum time problem on S 3 (14) > singular extremals are integral curves 
of the drift, i.e. they must correspond to a control almost everywhere vanishing. 

Since for a fixed u 6 [—1, 1] every trajectory of (14) is periòdic with period . 47r we have that: 

yu 2 sin 2 q+cos 2 a 

Proposition 11 Given an extremal trajectory 7 of type B t (resp St), then t < Air (resp. t < —^— ). 

The following proposition describes the switching behavior of abnormal and bang-bang normal extremals (see 
Section A.l for the definition). 

Proposition 12 Let 7 be an abnormal extremal of (14)- Then it is bang-bang and the time duration between 
two consecutive switchings is always equal to tt. In other words, 7 is of kind B s B 7r ...B 7r B t with s,t < tt. 
On the other hand, if 7 is a bang-bang normal extremal, then the time duration T along an interior bang arc 
is the same for all interior bang arcs and verifies tt < T < 2ir (i.e. 7 is of kind B s Bq-...Bq-B t with s,t < T). 

For the optimal trajectòries containing a singular arc we have the following: 

Proposition 13 Let 7 be a time optimal trajectory containing a singular arc. Then 7 is of the type B t S s B t ' , 
with s < 2v if t > or t' > and s < 47r otherwise. 

— cos a J cos a 

Thcsc results on S 3 ~ SU (2) are useful to determine the optimal synthesis on Sb , since every optimal trajectory 
011 is the projection of an optimal trajectory on S 3 . This is a simple consequence of the fact that Sb is an 
homogeneous space of SU(2): 
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Proposition 14 A time optimal trajectory 7 for the system (10)-(13) on Sb starting at Pjsr is the projection 
of a time optimal trajectory of (14) starting from a point satisfying \ipi\ 2 = 1 (recali that ip = (ipi,ip2) T G S 3 C 
C 2 ). 

Rcmark 15 Noticc that, since two opposite points on S 3 project on the same point 011 Sb ; it is easy to see 
from Proposition 11, that the projection on Sb of an optimal trajectory of (14) of type Bt (resp St), must 
be such that £ < 2tt (resp. t < -^ E — ). More precisely, for a fixcd u G [—1, 1] cvery trajectory of (10)-(13) is 
periòdic with pcriod 271 (the period divides by two after projection). 

1/ u 2 sin 2 q+cos 2 a 

B.1.2 Construction of the Synthesis on Sb 

In this section we assume et > 7r/4. Following Appendix A we first need to determine the sets A^ 1 (0), A^ 1 (0), 
and the function fg. Checking where Fs is parallel to Gs and where Gs is parallel to [F$,Gs], one gets 
= iv e &b '■ 2/2 = 0} and A^ 1 (0) = {y G Sb ■ V3 = 0}. To find the function fs we can choose 
for instance the coordinate chart defined on each hemisphere by the projection on the plain {(2/1,2/2) G K. 2 } , 
obtaining fs = (sin (1)2/3/2/2. Thcn Lemma 3 says that, every optimal trajectory bclonging to one of the regions 
{2/ G Ss : 2/3 > 0, 2/2 > 0}, {2/ G Ss : 2/3 < 0, 2/2 < 0} is bang-bang with at most one switching. Moreover only the 
switching from control —1 to control +1 is allowed. On the contrary, on the regions {y G Sb : 2/3 > 0, 2/2 < 0}, 
{y G Sb : 2/3 < 0, 2/2 > 0}, the control can switch only from +1 to — 1. Moreover, thanks to Lemma 1, every 
singular extremal must lie on the equator. The following lemma characterizes the structure of the bang-bang 
extremals for the problem (P). 

Lemma 5 Recali that t\ = ir — arccos(cot 2 a) and Í3 = 7r+arccos(cot 2 a) and consider a bang-bang extremal for 
the problem (P). Then it is of the form B S B V ^B V ^ . . . with s G [0, £1] U [ir, Í3], where, on the set [0,íi[U[7r, £a[, 
V (.) is defined as follows: 

. . / sins \ 

v(s) := tt + 2 arctan 5 — . 

\coss + cot z q / 

If a = 7r/4 then t\ = £3 = 7r and v(tt) := tt, while if a> tt/A we set v(t\) := «(£3) := 2tt. 

Notice that the function v(.) has the same expression (17) obtained in the case a < tt/A (excepted at the points 
£1 and £3). However its interval of definition is different. 

Proof of Lemma 5. As shown above, the meridian A^ 1 (0) and the equator A^ 1 (0) divide the sphcrc in 
four parts and in each of them the sign of the function fs is constant and changes when passing through 
A^ 1 (0) or A^ 1 (0). In particular, following j + or 7" (cfr Remark 4) in the case in which a > tt/A this happens 
at the times £1 (where the equator is crossed), at time tt (where A^ 1 (0) is crossed) and at time £3 (again is 
the equator to be crossed). Applying Lemma 3, we obtain that for an extremal trajectory the first switching 
may oceur only on the intervals [0, £1] and [71", £3]. Exactly as in [11], one shows that the extremal must have 
the form B s B v r s \B v r s \ . . . with s G [0, £1] U [tt, £ 3 ]. The case a = tt/A is similar. H 

Rcmark 16 One can also show that every trajectory starting from P/v, of the form B s B v r a \B v r s \ . . . with 
s G [0,£i] U [71", £3] is extremal i.e., for every s in such set, there exists an initial value of the covector A such 
that the switching function <p{.) vanishes for the first time at time s. 

Unlikc the case in which et < tt/A, in the case et > tt/A it is possible to establish the presence of optimal 
trajectòries containing a singular arc, whose switching strategy is described by the following proposition, 
illustrated in Figure 10 A. 

Proposition 15 Let a > tt/A. A trajectory 7 of (10)-(13) starting with control u = 1 and containing a 
singular arc is a solution of (P) if and only if it is of the form B t S s B t ' and satisfies the following conditions: 

• £ = £1 = tt — arccos(cot 2 a) i.e. 7 coincides with 7 + until it reaches the equator. 

• s < arccos(cot et)/ cosa i.e. the singular arc is optimal until it reaches the point + = (1,0, 0) r . 
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Figurc 10: The region covcrcd by optimal trajectòries with singular arcs and the (Y, K) 3 frame point 



• If s = arccos(cota)/cosa, then the trajectory is of type B t S s (i.e. the time duration of the last bang arc 
reduces to zero). If s < arccos(cota)/ cosa, then 7 is optimal until the last bang arc reaches the equator 
(i.e. it does not exist t s]0,í'[ such that 7(4 + s + t) is contained in the equator). 

An analogous result holds for trajectòries starting with control — 1. 

Rcmark 17 Notice that in the case a = tt/4, Proposition 15 provides a singular trajectory degenerated to a 
point. In other words for a = n/4 there are no singular trajectòries that are optimal. 

Rcmark 18 Notice that the previous result completely characterizes the optimal synthesis in some neighbor- 
hoods of the points = (±1,0, 0) T , namely U ílf, and moreover it determines the presence of two 
symmetric overlap curves contained inside the equator. The synthesis around the point + is represented in 
Figurc 10 A. 

Proof of Proposition 15. Consider a trajectory, that is a solution of (P) starting with u = +1 and containing 
a singular arc. Using Propositions 13 and 14 this trajectory must be of the form BtS s Bt> and, since the singular 
arc is contained inside the equator, we have t = t\ (the case t = Í3 can be easily excluded). Consider a 
singular arc containing in its interior the point + . This arc contains two points of the form (j/J, — y®, 0) T 
and (yJ,í/2iO) T j with both ï/J, y® positivc, that can be connected by a bang arc. Using classical comparison 
theorems for second order ODEs, one can easily compare the time needed to follow such trajectory with the 
time needed to steer the two points along the singular arc finding that the bang arc is quicker than the singular 
arc. Therefore a singular arc containing + cannot be optimal. By symmctry, the cxtremal trajectòries that 
have the same singular arc, but the last bang arc corresponding to opposite control, must meet on a point of 
the equator. Therefore the arc of the equator which is comprised between the point + (resp 0~) and the 
second intersection point with 7 + (resp. 7") is an overlap curve. It remains now to verify that the trajectòries 
described above are optimal (until the last bang arc reaches the equator). This is a straightforward consequence 
of the fact that the quickest bang-bang trajectòries that enter the region spanned by such trajectòries (i.e. the 
closure of the regions are not extremal because of Lemma 3. 

Rcmark 19 Notice the trivial fact that, if a trajectory 7 defined on the interval [a, b] is optimal between 7(0) 
and 7(6), then the restriction of 7 in [c, d], c, d £ [a, b], c < d, is optimal between 7(0) and j(d). 

Using Remark 19, we have that Proposition 15 characterizes completely the time optimal synthesis on 
P N A ± and in the closure of flf U flf, i.e. it proves items T1-T6, T9 and T10, of Theorem 1. 

Rcmark 20 From Lemma 5 we obtain that there are four familics of bang-bang trajectòries. In particular the 
famílies starting with control +1 and switching rcspcctively in [0, íi] and [n, Í3] join at the point B + , generating 
an amazing (Y, K)z frame point, in the framework of the classification of [12]. See Figure 10 B. 

Next we give the proof of Proposition 2, from which it follows Tll of Theorem 1, and, using again Remark 
19, also T7. 
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Proof of Proposition 2 By Proposition 15, there are no optimal trajectòries containing a singular arc 
joining Pm with Pg. One can easily see that the only possible trajectòries steering P/v to Ps with only one 
switching are those described in the statement of the proposition, that we have to compare with trajectòries 
having more than one switching. Trajectòries having two switchings with the first or the last bang longer than 
7r and trajectòries with more than two switchings arc exeluded since from Lemma 5 their total time is larger 
than 2ir. Trajectòries having two switchings and length of the first arc Si and the length of the last arc Sf 
satisfying Si,Sf < n are exeluded since by Lemma 4 they must satisfy Sj = Sf. For these trajectòries the total 

time can be easily computed and it is 2tt + 2 arcsin f a^^j j > 2tt. I 
ítem T8 is proved by the following: 

Proposition 16 If y € U {DPg \ Ps), then *Yg is made by a unique trajectory of the form B t B t ' , with 
< t < t\ and the first bang corresponding to control +1. A similar result holds if y € íl^ U {DPs \ Ps)- As a 
consequence there is not a cut locus in the region il^ U íl^ . On the other hand ^^asty ^ ^nasty contains a cut 
locus. 

Proof of Proposition 16 Dcfine the following three families of extremal trajectòries: 

yf(t) := e tx s e sX sP N , with s e]0,h[ and t<v(s), 
7f(t) := e tx s e sX sP Nl with s € [ir,t 3 [ and t<v(s), 
7 f(í) : = e tx s e v{s)x se sX s P Nl with sG]0,íi[ and t<v(s). 

First notice that from Proposition 2, there are no optimal trajectòries of kind 7^ reaching the arc BPg . 
Now for every point x <E DP$ the following happens: i) there exist sa, t>A such that x — 7^,(ía)j and they 
are unique; ii) if there exist Sb, ïb (resp. sc, te) such that x = 7^(íb), (resp. x = 7^(íc))j then they 
are unique. By direct computation, one can compare the times the three trajectòries need to reach x, i.e. 
sa + ía, sb + íb, sc + v(sc) + tc, finding that the optimal trajectory is of kind j A (these computations are 
long, not very instructive, and we omit them). From this fact the first part of the claim immediately follows. 
Morcovcr it implies that there is not a cut locus in Í2^~, since the only trajectòries entering such region are 
those of the form r y A . The existence of a cut locus in f^asty ^ s evident, since no optimal trajectory bclonging 
to the families 7^, 7 B , 7° leaves íínasty ^ ne re asoning in and in íín asty is similar. H 



End of the proof of Theorem 1 

To conclude the proof of Theorem 1. it remains to prove T12. Considcr by contradiction an optimal bang- 
bang trajectory 7 defined in [0,í 7 ] steering Pjv to a point of í^asty w ith at lcast three switchings. Define 
t = max{í G [0,í 7 ] : 7(í) ^ ^nastyi- Then, by Remark 19, j\[o,í\ must be optimal between P/v and 7(í). Then, 
from the results proved above, we deduce that 7|[q ï] can have at most one switching. Therefore 7 switches at 
lcast two times in íí^asty an< ^ ^ nc arc between them must be complctely contained in í^n ast y and this lcads to 
a contradiction since the sign of fs is constant in íí^astv ( see Lemma 3). H 



Bcfore proving Proposition 3, notice that the point D + , which is obtained following the trajectory 7+ for 
a time 7r (sec Figuro 6), belongs to two diffcrent families of bang-bang trajectòries at time w, one given by 
trajectòries starting with control —1 and switching at time s < íi, the other one given by trajectòries that start 
with control 1 and switching at time s € [ït, Í3]. Moreover, since v(0) = 7r, there must be a switching curve 
starting at D + and therefore we deduce that there are two possible behaviors of the optimal synthesis around 
this point: either this switching curve is optimal or the two fronts continue to intersect generating an overlap 
curve. 

Observe that if a > 7r/3 the trajectòries of the type B 8 B v r s \Bt with s small cannot be optimal since the vector 
fields Xg and Xg point to oppositc sides on the switching curve (i.e. the switching curve "refleets the trajec- 
tòries", see footnote 1). In this case the two families of bang-bang trajectòries described above must intersect 
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giving rise to an ovcrlap curve. Thereforc to provc Proposition 3 wc assumc a < n/3. 

Proof of Proposition 3 First we parameterize the switching curve with respect to the first switching time 
(assuming without loss of generality that this curve starts with u = —1): 

C(s) = e v(s)x s e sX s P N . 

We consider the functions £i(s) = det (C(s), C'(s), Xg (C(s))) (here the superscript ' denotes the derivative 
with respect to s) and £2(5) = det (C(s),C'(s),Xg(C(s))) . It is easy to see that the optimality of C(.), for s 
small, depends on the signs of such functions. Indeed C(.) is locally optimal near the point D + = C(0) if and 
only if for evcry s > (small cnough) and given a neighborhood of C(s) which is divided in two connected 
components Ui, U2 by the trajectory C(.), both Xg(C(s)) and Xg(C(s)) point towards Ui or towards C/2. It 
is easy to see that this happcns if £i(s) and £2(5) have the same sign. Notice that £i(0) = £a(0) = and that 
£i(s) = det (P/v, Xg{Pif), e~ sX s X$ (e sX s Pjy)) = 2cosasin 2 asins, which is positive for every q < 7r/2 and 
s €]0, 7r[. To determine the sign of ^(s) near it is enough to look at the sign of the derivative Ç' 2 (0) which can 
be computed directly: Ç' 2 (Q) = 4cosasin 2 a(l — 2sin 4 a). We deduce that, if a < arcsin(l/v / 2), the switching 
curve C(.) is optimal for s small enough. For the particular value a = arcsin( 1/^/2) one can easily check that 
the function ^(0 is negative for s > small, and then C(.) is no more optimal for a > arcsin(l/^2). The 
tangeney of the switching curve starting at D + if a > n/A, is a consequence of the fact that, in this case, the 
bang-bang trajectory switching at D + is an abnormal cxtrcmal (sec Proposition 2 and [12], Proposition 23 pag. 
177). ■ 



B.2 Time optimal trajectòries reaching the south pole for at < ir/4 

The purpose of this section is to characterize the optimal trajectòries steering P/v to Pg in the case a < 7r/4, 
i.e. to prové Proposition 4 and 5. A key tool is Lemma 4. Recali the shape of the function v(s), in the case 
q < n/4 (see Figure 7 A). Given a < n/4 and s £ [0, n] with s ^ arccos(— tan 2 a), there exists one and only one 
time s*(s) € [0,7r] different from s, such that v(s) = v(s*(s)). From Section 3.2 recali the following definition 
of candidate optimal trajectòries: 

• Sf = s*(sí) (i.e. TYPE-l-candidate optimal trajectòries), 

• st = Si (i.e. TYPE-2-candidate optimal trajectòries) 

A uscful rclation bctwcen s and s*(s) is given by the following: 

Lemma 6 For a < n/4 and s £ [0, n], it holds s + s*(s) = v(s). 

Proof of Lemma 6 Both s and s*(s) satisfy the following equation in t £ [0, n]: 

cot Q U(S )) = - cos(í) Sl ; ( c ^ 2(a) =► cos Q,( S ) - t) = - cos Q,(.s)) cot 2 (a). 

Therefore, since ^v(s) — t £ [— n,n] Vs,í £ [0,n] and s*(s) 7^ s, it must be: s*(s) — h v ( s ) = 2" w ( s ) — s ^ 
s + s*(s) =v(s). ■ 

The description of candidate optimal trajectòries is simplified by the following lemma, of which we skip the 
proof. 

Lemma 7 Set: 

, / cot (|u(s)) — sin(a) 

Z(s) = - -cot (|«(s)) 

p \ sin(a) 

where p = ^ cot 2 (\v{s)) +sin 2 (a) . Then, if 6(s) is defined as in (20), we have e e ^ z ^ = e v ^ x s e v ^ x s . 
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Notice that the matrix Z(s) £ so(3) is normalizcd in such a way that the map 1 1— > e' 2 ^ £ 50(3) represents 
a rotation around the axes iï(s) = (0, sin(a), cot(^u(s))) with angular velocity equal to one. 

To prové the results stated in Section 3.2 we study separately the two possible cases listed above: 

Proof of Proposition 4. In this case we consider TYPE-l-candidates optimal trajectòries. Assume that 
the optimal trajectory starts with u = — 1 (the case u = 1 is symmetric) and has an evcn number n of 
switchings. Then it must be 

Ps = e SfX * e v(s * )x t. . . e v( - s ^ x ^ e s ' x s P N (29) 

n—X times 

where P/v and Ps denote respectively the north and the south pole, and we have that 

e SiX~ p s = e v( Si )X~ e v( Si )X+ _ _ _ _ e v( Si )X+ e SiX- p N = e ínO( Si )Z( Si ) eSi Xs 

from which we deduce that sí must satisfy 

—n9(si) = 7r + 2pn for some integer p. 

It is easy to see that a value of Sj which satisfies previous cquation with p > doesn't give rise to a candidatc 
optimal trajectory since the corresponding number of switchings is larger than Nm- Thercforc in previous 
cquation it must bc p = 0. If n is odd the relation (29) becomes 

Ps = e *f x t e v[St)x A.. . . e v{s * )x s, e s ' x s P N (30) 

tl— 1 times 

and, moreover, by symmetry: 

n = e 1 s e y ' . •?.... e v " s e s Ps- 
Then, combining with (30) and using the relation Lemma 6, we find: 

P N = e ~ SiX s e v( ~ Si)x .s. . . . e vi -^ )x s e^ x s P N = e -^ x s e ne{s * )z{s ^ e s ^ x s P N . 

2n times 

Since e SiX s P N is orthogonal to the rotation axis R(sí) corresponding to Z(sí), previous identity is satisfied if 
and only if n9(si) = 2mir with m positive integer. As in the previous case, for a candidate optimal trajectory, 
it must be m = 1. 



Proof of Proposition 5. Hcrc wc consider TYPE-2-candidate optimal trajectòries. For simplicity call 
Si = Sf = s. Assume, as before, that the optimal trajectory starts with u = —1 . If this trajectory has 
n — 2q + 1 switchings then it must be 

P s = e sX s e qe{s)z{s) e sX s P N . 

In particular the points e~ sX s P s and e sX s P N must belong to a plane invariant with respect to rotations 
generated by Z(s) and therefore the difference e sX s P N — e~ sX s Ps must be orthogonal to the rotation axis 
R(s). Actually it is easy to see that this is true for every value s £ [0, n], since both e~ sX s Ps and e sX s P N arc 
orthogonal to R(s). Moreover, since the integral curve of Z(s) passing through e sX s P^ and e~ s P s is a circle 
of radius 1, it is easy to compute the angle /3(s) between these points. In particular the distance between e sX s P N 
and e~ sX s P s coincides with 2 sin(^p-) , and so one easily gets the expression /3(s) = 2 arccos(sin(a) cos(a)(l — 
cos(s))). Then Proposition 5 is proved when n is odd. 

Assume now that the optimal trajectory has n = 2q + 2 switchings, then we can assume without loss of 
generality that P s = e sX s e v ^ x s e q9 ^ z ^e sX s P N . First of all it is possible to see that e~ v ^ x Í e~ sX s P s 
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is orthogonal to R(s). So it rcmains to computo the angle (3(s) bctween the point e sX s P/v and the point 
e -v(s)x+ e -sX s p s Qn ^ c plane orthogonal to R(s). As before the distance between these points coincides 

with 2sin(^p-). Instead of computing directly /3(s) we compute the difference between the angles /3(s) and the 
angle /3(s). We know that 

2sin(^-/3(s)) = \e- v ^ x s e~ sX s P s - e - sX s P s \ = \e~ sX s P s ~e v ^ s)x s e~ sX s P s \ = \e~ sX s P s ~ e s '^ x s P s \. 
Using the fact that s and s*(s) satisfy the relation v(s) = v(s*(s)) one can easily find that 



\ e - sX sP s - e s *^ x sP s \ = 2^1 - cos 2 (a)sin 2 

Therefore /3(s) = $(s) + 2 arecos (cos(a) sin (|u(s))j . This leads to f3(s) — (3(s) = 0(s)/2 and the proposition is 
proved also in the case n is even. H 



B.2.1 Proof of Proposition 6, on the alternating behaviour of the optimal synthesis 

In this section we need to consider also the dependence on a of the functions v(s), 8(s), (3(s), ^(s), G(s). There- 
fore we switch to the notation u(s, a), 8(s, a), /3(s, a), T{s, a),Ç{s,a). 

The claims on existence of solutions of Proposition 6 come from the fact that J-(Q) = T{tt) = and the 
only minimum point of T oceurs at s = tt — arccos(tan 2 (a)). It turns out that the image of T is a small interval 
whosc lcngth is of order a and therefore equation (22) has a solution only if a is close enough to ^ for some 
integer number m. This proves C. i.c. the existence of ^(m) satisfying ^(m) = 0(l/m). 

On the other hand it is possible to estimate the derivative of Q with respect to s showing that it is negative 
in the open interval ]0, 7r[. Therefore, since G{0) = ^ -I- 1 and G{ir) = — 1, equation (23) has always two 
positive solutions. 

For the particular vàlues a = where m > 1 is an integer number, the solutions to the equations (22) 
and (23) give rise to two candidate optimal trajectòries, one with m bang arcs, all of length tt (TYPE-1 and 
TYPE-2 candidate optimal trajectory at the same time), while the second one has one more switching and 
is a TYPE-2 candidate optimal trajectory. We want to see that the optimal trajectory is the first one. To 
this purpose, we need to estimate the time needed to reach the south pole by the second candidate optimal 

2 

trajectory showing that it is greater than mir = ^j. 

First, using the Taylor expansions with respect to at and centered at of /3(vr/2, a) and 6(7r/2,a), one 
obtains 

= ^-«7 + «(«)■ (31) 
2 2a 4 

We want now to estimate the solution s(a) of the equation G{s, a) = This can be done using the previous 
expression (31) and since it is possible to estimate the derivative of G(-) with respect to s near s = n/2, in the 
following way: 

±g( a ,a) = -l + o(\a\ + \^ 
thcn it is casy to find that s(a) — — — a— + o(a), and, consequently v(s(a), a) — n + 2a 2 + o{a 2 ). Therefore 



- - s 



2s{a) + í 1 ) v(s(a), a) = — h a— + o(a). In particular, for a = ^ this expression coincides with the 



f TT \ TT" TT 

(__lj v(s(a))a) = _ +Q; _ 

time needed to reach the south pole by the candidate optimal trajectory and, since for m large enough it is 

2 

larger than rmr = |— , we found that this trajectory cannot be optimal. Since the solutions to the equations 
(22), (23) change continuously with respect to a for each fixed number of switchings n, wc easily deduce that 
if wc slightly deercase a starting from the value the solution of (22) for n = m does not give rise to an 
optimal trajectory. 

For a slightly smaller than a :— there is a TYPE-2 candidate optimal trajectory corresponding to a 
solution (si(a),m + 1) of (23), where Si(.) is continuous (on [a — s,a\) and si(a) = 0, and there is also a 
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TYPE-1 candidate optimal trajectory corresponding to a solution (s2(a),m) of (22) where S2(.) is continuous 
(on [ce — e,a]) and s 2 {a) = 0. Clearly for ce = ce these trajectòries coincide. So we have to compare the time 
to reach the south pole for such trajectòries with ce close to a. 

We start with the TYPE-1 candidate optimal trajectory. From equation (22) we have that —Ofaice), ce) = 0. 

da 

We use a subscript s, a to denote the partial differcntiation with respect to such variables. Since S (O, ce) = 
we cannot apply directly the implicit function theorem near (0, ce). However, if we set §2(0) = s 2 (ce) we find 



that s' 2 {ce) 



2s 2 (a)9 a (s 2 {a),a) 



(the supcrscript ' denotes differcntiation with respect to ce), and then, passing 

2 



6 s (s 2 {a),a) 

to the limit as (s2(oe),a) tends to (0,a), one easily finds that s'nia, - 

sn^ay cos(qí) 

Now we want to determine the way in which the total time T2(a) = mv{s2{ce), ce) changes. It is easy 
to see that T2(ce) is not diffcrcntiablc at ce, therefore we introduce the function F(a) = (T 2 (a) — T2(ce)) 2 = 
m 2 (v(s2(ce) 7 ce) — 7r) 2 . 

Then — F(a) = 2m 2 — v(s2{oe), ce) (i>(s2(a), ce) — n) = 2m 2 (v s (s2(ce), a)s' 2 {ce) + v a (s2(ce), ce)j and, after 



the substitution s' 2 (ce) = 



2s 2 (a) 



we can pass to the limit as ce converges to a obtaining 



da 



F(a) 



m 2 v 2 (0, a)s' 2 (a) = — 8m 2 tana. 



Now we consider the TYPE-2 candidate optimal trajectory and we want to estimate si(a). From equation 
(23) we have that si(.) is implicitly defined by the equation <ï>(si(a),a) := 2/3(si(a), ce) — m6(si(ce),ce) = 0. 
As before it is easy to see that s±(.) is not differentiable at ce and therefore we introduce the parameter 
Si(ce) = s 2 (ce). As before, it is possible to compute the derivative s'i{ce): 



Sí (a) 



lim 

a — >a 



2s 1 (ce)<è a (s 1 (a),ce) 



2m 



$ s (si(a),a) sinacosa(l + msin a) 

We have now to estimate the total time ï\(a) = 2s2Íce) + mz)(s2(a), ce) for ce close to ce. Dcfine 

G(ce)= (T^ai-T^a)) 2 = (2s 2 (a) + m(v(s 2 (oe), a) - tt)) , 
then we compute the derivative of G(.) as follows: 



Ta G{a)\ a , 



lim 

OL — >Oi 



lim 

a — >q 



2 2 + m 



^2s 2 (a) + m(v(s2(ce), ce) - tt) 
w(s 2 (a), ce) — v{0, a) 



82(0) 



v a (s2(ce),ce)s 2 (a) 



+ v a (s 2 (ce),ce) 



S2(ce) 



lim 

OL — >Oí 



= (2 + mv s (0, ce)) s' 2 (ce) = -(2 + 2msin 2 a) 



2s 2 (a) 

s' 2 (ce) + m(^v s (s 2 (ce) 7 ce)s 2 (ce) + v a (s 2 (a),a)s 2 (a) 

8m(l + msin 2 oi) 



2m 



Since 



sin ce cos a(l + m sin 2 ce) 
8m(l + m sin 2 ce) 



sm ce cos ce 



> mtana 



sm ce cos ce 

we deduce that G(a) decreases faster than F(ce) as ce goes to à and, since Xi(a) and T 2 (a) are decreasing for 
ce close to ce, we have that T 2 (a) > T\{ce), i.e. the TYPE-1 trajectory is optimal for ce € [ce — e,ce\. 
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C The Time Needed to Reach Every Point of the Bloch Sphere 
Starting from the North Pole in the case a 6 [7r/4, 7r/2[ 

In this section we assume a 6 [7r/4,7r/2[. If a is close to 7r/4 it is easy to verify that the south pole is not 
the last point reached by bang-bang trajectòries (the last point reached belongs to the cut locus present in the 
region í^nasty) an d the time needed to cover the whole sphere is slightly larger than 2ir. 

On the other hand, if a is large enough then the velocity along a singular arc is small and therefore the time 
needed to move along trajectòries containing singular arcs is larger than 2tt. The following proposition gives 
the asymptotic behaviour of the total time needed to reach every point from the north pole and determines 
the last point reached by the optimal synthesis for a large enough. 

Proposition 17 Let T(a) the time needed to cover the whole sphere. Then, if a is large enough 

„, , tt 2 arcsinfcot a) . tt „/7r \ , 

T(a) = h tt + 2 arcsinícot 2 a) = h tt - 2 + O a) (32) 

2 cos a cos a 2 cos a V 2 / 



and the last points reached for a fixed value of a are ±(vi — cot 2 a, cot a, 0) T . 

Proof of Proposition 17 From Proposition 2 the last points reached by optimal trajectòries of the form 
B t S s B t i must lie on overlap curves which are subsets of the equator. Therefore it is enough to estimate the 
maximum time to reach these overlap curves. Assume that the frrst bang arc corresponds to the control u = 1 
and denoto by (3 the angle corresponding to the arc of the equator between the last point of the singular arc 
and the point + = (1,0, 0) T . Notice that j3 s]0, arccos(cota)[. Then it is easy to find the expression T(a,/3) 
of the time needed to reach the overlap curve along that optimal trajcctory: 

o arccos(cot a) (3 ( cos 2 a — tan 2 (3 

T(oí, p) = tt — arccos(cot a) H h arccos 1 



cos a cos a \ cos 2 a + tan (3 

The conclusion follows finding the maximum with respect to [3 of the previous quantity, which corresponds to 
the value (3 = arcsin(cot ce). Notice that (3 belongs to the interval of defmition of (3 only if a > arccot(\/2/2). 



Remark 21 If ot > arccot(v / 2/2) there are two symmetric neighborhoods of the points ±(-\/l — cot 2 a, cot a, 0) T 
that are not reached by optimal trajectòries at time t < T{a) with t close enough to T(a), i.e. the reachablc 
set is not simply connected. 

Remark 22 Recali that for system (6) the time needed to cover the whole sphere for a close enough to tt/2 
is obtained dividing by k = c 2 ^ a the expression (32). Therefore, if we fix E it turns out that this quantity 
converges to ^ as M goes to infinity. 
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